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SUMMARY 

A nondimensional model of microstructural ly short crack growth in creep 
is developed based on a detailed observation of the creep fracture process of 
304 stainless steel. In order to deal with the scatter of small crack growth 
rate data caused by microstructural inhomogeneity, a random variable technique 
is used in the model. A cumulative probability of the crack length at an arbi- 
tary time,_G(a,t) , and that of the time when a crack reaches an arbitary 
length, F(t,a), are obtained numerically by means of a Monte Carlo method. 
G(a,t) and F(i,a) are the probabilities for a single crack. However, multi- 
ple cracks generally initiate on the surface of a smooth specimen from the 
early stage of creep life to the final stage. Taking into account the multi- 
ple crack initiations, the actual crack length distribution observed on the 
surface of a specimen is predicted by the combination of probabilities for a 
single crack. The prediction shows' a fairly good agreement with the experimen- 
tal result for creep of 304 stainless steel at 923 K. The probability of creep 
life is obtained from an assumption that creep fracture takes place when the 
longest crack reaches a critical length. The observed and predicted scatter of 
the life is fairly small for the specimens tested. 


INTRODUCTION 

In designing mechanical components for high temperature applications, the 
engineer must consider the failure mechanism at such temperature. Accumula- 
tion of creep damage often causes failure at high temperature. The creep 
damage usually consists of multiple cracks or cavities. For creep life predic- 
tion, it is necessary to evaluate the damage accumulation based on the synthe- 
sis of crack growth and/or cavity growth. 

Creep fracture mechanisms were experimentally investigated in detail using 
304 stainless steel smooth specimen in previous works (refs. 1 to 4). Observa- 
tions revealed that: (1) cracks continuously initiated at the surface from the 

early stage of life to the final stage; (2) the cracks appear to grow randomly 
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randomly owing to the effect of local mi crostructural inhomogeneities; and (3) 
crack coalescence and necking of the specimen (reduction of area) accelerated 
drastically the growth in the final stage of life. Eventually, the creep life 
is governed by the crack initiation and growth. 

The existing fracture mechanics law for long creep crack growth (refs. 5 
to 11) is not directly applicable to the random growth of mi crostructural ly 
short cracks in creep. Few studies have been reported on the short creep crack 
growth (ref. 12) because of the difficulties in interpreting the randomness. 

The authors previously proposed a model of the mi crostructural ly short creep 
crack growth (ref. 12). In the model, a random variable technique was used to 
deal with the distribution of crack growth data. 

In this paper, the model is recast into a nondimensional form for predic- 
tion of the creep crack growth. This greatly reduces the numerical calcula- 
tions required. Taking into account the continuous initiation of cracks and 
the failure crack length the cumulative probability of creep life is predicted 
by a Monte Carlo simulation method on the basis of the nondimensional model. 


NOMENCLATURE 


a half crack length 

a c half crack length at which the fluctuation of crack growth rate 

blends into long crack behavior (fig. 1(b)) 

af critical half crack length for creep rupture 

B coefficient of creep power law 

C coefficient of the relation between crack growth rate and crack 

length in the long crack growth law 

C c coefficient of the relation between crack growth rate and creep 

J-integral in the long crack law 

D grain boundary length between two adjacent triple points 

D av average of D 

da/dt crack growth rate 

( da/d t ) |_ crack growth rate at a triple point (fig. 1(b)) 

(da/dt)(j crack growth rate at a medium point between adjacent triple points 

(fig. 1(b)) 

F(t,a) cumulative probability of f(t,a) (probability that a crack reaches 
length a before t) 

t 

f(t,a) probability density function of a time t at the instant a crack 
reaches length a 
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G(a, t) 

g(a,f) 

h 

J* 

K(t) 

k(t) 

Ll 

l 2 

Mj 

m 

mi 

m2 

n 

P(t r ) 

S 

t 

ti 

tr 

a(n) 

|3 and y 

At 

At] 

a 

a g 


cumulative probability of g(a,t) (probability that crack length is 
shorter than a at t) 

probability density function of a crack length a at time t 

nondimensional length between crack tip and grain boundary triple 
point (fig. 2) 

creep J-integral 

number of cracks per unit surface area (crack density) 
time rate of change of K(t) 

constants in the relation between crack growth rate and crack length 
for a small crack (fig. 2) 

constants in the relation between crack growth rate and crack length 
for a smal 1 crack (fig. 2) 

boundary - crack shape correction factor 

number of cracks at t = t r (= K(t r )S) 

constants in the relation between crack growth rate and crack length 
for a small crack (fig. 2) 

constants in the relation between crack growth rate and crack length 
for a small crack (fig. 2) 

exponent of power creep law 

cumulative probability of creep life 

specimen surface area 

time 

crack initiation time 
creep life 

function of n (Eq. (4)) 

positive integers 

l/k(t r )S (Eq. (17) or Eq. (A-3)) 

t r /(J (Eq. (A-4) ) 

stress 

applied stress 
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<Ka,t) cumulative probability of a crack length a for multiple cracks case 
at a time t 

Note: bar over variables indiates nondimensional form. 


SMALL CRACK GROWTH MODEL 

Figures 1(a) and (b) show schematic growth curves and growth rate of 
microstructural ly short cracks initiated on the surface of a smooth specimen 
of 304 stainless steel under creep conditions (ref. 3). The characteristics 
of the crack growth are summarized as follows. 

(1) Cracks initiate at grain boundary triple points or carbide precipi- 
tates and grow at high rates along grain boundaries but decelerate as they 
approach triple points or sharp bends of grain boundaries. 

(2) The relation between crack growth rate and crack length is greatly 
dependent on individual cracks. Fluctuation in crack growth rate is large for 
cracks smaller than a few multiples of the grain size, a c . 

(3) The degree of fluctuation diminishes as the crack length becomes 
longer and blends into the scatter band for "long-crack" behavior. 

A stochastic model of microstructural ly short crack growth was proposed 
based on the experimental results in a previous paper (ref. 12). The model is 
summarized briefly as follows. 

(1) The crack length, a, and the grain boundary length between two adja- 
cent triple points, D, are the projected lengths on the plane perpendicular to 
the applied stress axis. 

(2) The grain boundary length between two adjacent triple points, D, is a 
random variable with a normal distribution. 

(3) The intergranular crack growth rate alternates between an upper and a 
lower peak depending on the location of crack tip as shown in figures Kb) and 
2 . 


(4) The upper and lower peaks of the crack growth rate are modeled by ran- 
dom variables of logarithmic normal distributions. Their medians and 90 per- 
cent confidence bounds are listed in table I. This procedure gives reasonable 
crack growth data as compared with the experimental curves (ref. 12). 

(5) The crack growth rate for a long crack is proportional to the crack 
length. 


da 

dt 


= Ca . 


( 1 ) 


(6) The crack growth rate between grain boundary triple points is given by 
the equations shown in figure 2. 
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The creep crack growth rate, da/dt, has a good correlation with a creep 
.J-integral (C*-parameter or modified J-integral), J* (refs. 5 to 11). For many 
materials, da/dt is proportional to J*. 


da _ 
dt “ 



( 2 ) 


where C c is a material constant (ref. 13). J* for a finite crack in an 
infinite plate is evaluated by (refs. 1 and 14) 

J* = M,a(n)Bo" +1 a (3) 

J 9 

a(n) = 3.85 <n - — + (4) 

\/fi n 

where Mj is the boundary - crack shape correction factor, og is the applied 
stress, and B and n are the coefficient and the exponent of creep power law 
(Norton's law), respectively. Substituting equation (3) into equation (2), we 
obtai n 


4f = C M,a(n)Ba" +1 a 
dt c J g 

= Ca . (5) 

The crack length, a, the time, t, and the stress, o, are nondimensional- 
ized as follows, 


a 



( 6 ) 


t = C o Bo"t 

egg 


(7) 


o 





( 8 ) 


where D av is the average of the grain boundary length between two adjacent 
triple points. On the basis of equations (6) to (8), the relation between the 
nondimensional crack growth rate, da/dt, and the nondimensional creep 
J-i ntegral , J* , gives 


Ca (9) 

dt 


where 


C = Mj<x(n) 


( 10 ) 
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It should be noted that the nondimensional crack growth law is independent of 
the applied stress. The random variables used for the upper and lower crack 
growth rates are listed in table I while the normalized variables are listed 
in table II. We are also able to normalize the quantities, D, h, m^ , l_i , m 2 , 
and L 2 shown in figure 2. Therefore, any of the results can be derived from 
the nondimensional equations. Since theresults arei ndependent of the applied 
stress, both of the probabilities of G(a,t) and F(t,a) at any stress and at 
any time (to be described in the following section) can be expressed by the 
nondimensional analyses. Results can be converted back into the dimensional 
form through equations (6) to (8). 


PROBABILITY OF CRACK GROWTH 

The crack growth model is stochastic so that the growth curve has a proba- 
bilistic distribution as shown schematically in figure 3. g(a,t) is a proba- 
bility density function of crack length at a time t after crack initiation, 
and f(t,a) is aprobabi 1 i ty density function of time at the instant a crack 
reachesjength a. Their cumulative probabilities are represented by G(a,t) 
and F(t,a). G was partially discussed in the previous work (ref. 12), how- 
ever, F is introduced and discussed here for the first time. 

G(a,t) and F(t,a) are numerically given by the following procedure 
(Monte Carlo method) as shown in figure 4 on the basis of the model: 

(1) Generate random numbers. 

(2) Determine random variables D, (da/dt)y and (da/dt)|_. 

(3) Calculate the a - t relation for a crack. 

(4) Repeat 1 to 3 and get the a - t relations for many cracks. 

Figure 5 shows the change of G(a,f) calculated by a simulation of 300 
cracks. The constants used in the simulation are those derived from the creep 
behavior of 304 stainless steel at 923 K (refs. 2, 11, and 12) and are listed 
in table III. The convergence of simulation results wasconfirmed by a prelim- 
inary calculation except at the two extremes where G(a,t) is nearly equal to 
zero or one (ref. 12). G(a,t) can be physically interpreted as the probabil- 
ity that a crack length is shorter than a at time t. It becomes clear from 
figure 5 that the distributions change verylittle within the region G(a,t> < 
0.5, while they are strongly depending on t for G(a,t) > 0.8. This result 
suggests that the first grain boundary triple point has a high probability of 
blocking the crack growth and that the crack grows faster once it passes 
through the triple point. Similar tendencies were found in experiments 
(refs. 1 and 3). 

Figure 6 shows F(t,a) which is obtained by thesimulation of 30 000 
cracks in order to clarify the distribution near F(t,a) = 0. F(t,a) is the 
probability that the time when a crack length reaches a is shorter than t. 
The distribution shows that there is a difference in the shorter t region, 
while in the longer t region differences disappear. The difference at the 
shorter t' affects the probability of creep life. 
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CRACK LENGTH DISTRIBUTION 


Multiple cracks generally initiate continuously on the surface of a smooth 
specimen during the entire creep test. The actual distribution of crack 
lengths observed at a given time t is brought on not only by the distribu- 
tion of crack growth rates but also by the difference in the crack initiation 
time. G(a,t) obtained in the previous section is the cumulative probability 
of crack length at t after the initiation of a single crack (or that for mul- 
tiple cracks initiated at the same time). Therefore, it is necessary for the 
prediction of the actual crack length distribution to be coupled with the dis- 
tribution for crack initiation. 


The cumulative probability for an actual crack distribution, <Ka,t), at a 
time t is predicted by the following equation: 


4>(a,t) 


't 

G(a,t - t.) k(t.)dt./K(t) 

Jo iii 


( 11 ) 


where K(t) is the number of cracks per unit surface area (i.e., crack densi- 
ty) at time t, ti is a nondi^mensional^ crack initiation time, and k(t) is 
the time rate of change of K(t). <Ka,t) canbe interpreted as the probabil- 
ity that the crack lengths are shorter than a at a time a. The numer of 
cracks initiated from time t to t + dt is given as 

dK(t) = dt = k(t)dt . (12) 

dt 


Experimental observations have shown that K(t) was nearly proportional to t 
for conditions of creep (refs. 1 and 2). Thus, k(t) is constant and equation 
(11) can be written in the form 

t 

G(a,t - t. )dt. . (13) 

0 


4>(a,t) = 1 
t 


The lines shown in figure 7 are the predicted probabilities of crack 
length for a smooth specimen of 304 stainless steel at 923 K. Figures 7(a) 
and (b) are the probabilities at the medium creep life and near the end of 
creep life, respectively. The figures also include the actual distribution 
obtained from a creep test in air (ref. 2). The predictions agree well with 
the test results. Similar successful results are obtained at other, times t 
in air as well as the experimental results (ref. 2) in vacuum. 


CREEP LIFE OF SMOOTH SPECIMEN 

Creep tests of 304 .stainless steel revealed that the creep failure of 
smooth specimens took pl’ace soon after the longest crack reached a critical 
crack length, af (ref. 2). Crack coalescence and tertiary creep (reduction of 
area) brought on the drastic acceleration of crack growth and the unstable 
fracture. For the tests conducted in the previous work (ref. 2), the value 
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of af was about 0.15 mm. This implies that the growth of microstructural ly 
short cracks govern the majority of the creep fracture. 

Assuming a critical crack length for rupture, the cumulative probability 
of creep life, P(t r ), becomes equivalent to the probability that the longest 
crack reaches the length af. Here, multiple crack initiations must be taken 
into consideration. The derivation of P(t r ) from F(t,a)_is described as 
follows. Theprobabi 1 i ty that a crack initiated at time ti does not reach 
af at time t is given by [1 - F (t - tj , af)]. Then, 
m 

n {1 - F(t - tj.af)} is the probability that none of the cracks initiated at 
3 = 1 

1 1 , t£ , ... t m , reach af at t. PC t r ) , therefore, is given by 

m 

P(t r ) = 1 - n {l - F(t f - tj ,a f )| (14) 

3=1 

m = K(t r )S (15) 

Again, assuming k(t) is constant, equation (14) is converted into 


m 

p<y - ' - IT 1 - F ((m - j) • At ,a f ) ( (16) 

3=1 ’ 

At = (17) 

k(t r )S 


where S is the specimen surface area. 

P(t r ) is sensitive to the distribution near F(t,af) * 0, as shown in fig- 
ure 8 wherein the detailed distribution is given by the simulation of 150 000 
cracks. P(t r ) shown in figure 9 is obtained from the distribution of figure 8 
and equation (16). To simplify the solution of equation (16), an approximate 
estimation is used as described in the appendix. Figure 9 reveals that as the 
number of cracks, K(t r )S, increases, the creep life decreases and the scatter 
becomes smaller. In other words, the specimen surface area available for 
cracking will affect the creep life. Figure 10 shows the relationship between 
the nondimensional creep life and the number of cracks at P(t r ) =0.1, 0.5, 
and 0.9. The life changes considerably in the range of K(t r )S < 3000, but 
little in the range of 5000 < k(tr)S < 10 000. In the creep tests carried out 
previously, K(t r ) is about 20 mm -2 independent of the applied stress (ref. 2). 
Hence, the scale of S ; can be given and is shown on the upper side of 
figure 10. The real time scale at several stresses are also drawn at the 
right side of the same figure. The scatter of creep life is very small in the 
range of 250 mm 2 < S < 500 mm 2 . 
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S is equal to 320 mm 2 for the tested plate specimen of 304 stainless 
steel (width, thickness, and gauge length are 12, 4, and 10 mm, respectively) 
so that K(t r )S becomes 6400. The band in figure 11 shows the relationship 
between applied stress and predicted creep life for 0.1 < P(t r ) < 0.9 and 
K(t r )S = 6400. Also shown in figure 11 are the actual experimental results of 
creep rupture data (refs. 2 and 15). The predictions show good agreement with 
the data. 

The growth of surface cracks governs creep fracture in this relatively 
short time and high stress regime. However, it is expected that internal 
cracking plays an important role in the fracture process at lower stresses and 
longer times. Life prediction at lower stress may have to be modified to 
account for internal initiation and crack growth. 


CONCLUSION 

A nondimensional model of mi crostructual ly small crack for creep life pre- 
diction is proposed. The probability distributions of crack lengths and creep 
life are predicted based on the model. The results obtained are summarized as 
fol lows. 

1. The stochastic model of small crack growth is generalized. In the 
model, the scatter of crack growth rate owing to mi crostructual inhomogeneity 
is represented by means of random variables. The nondimensional model gives a 
general probability of crack growth independent of applied stress. 

2. The cumulative probability of crack length at a time t.GIa.t)^ and 
the cumulative probability of time when a crack reaches length a, F(t,a), are 
calculated by means of the Monte Carlo simulation method. 

t 

3. Taking into account the multiple crack initiations, the cumulative 
probability of crack length for the multiple cracks, 4>(a,t), is obtained from 
G(a,t). The predictions agree well with the actual distribution of crack 
lengths observed on a specimen of 304 stainless steel tested in creep at 923 K. 

4. The cumulative probability of creep life, P(t r ), is derived by assum- 
ing that the longest crack is responsible for the ultimate fracture. Predic- 
tions and experimental results are in good agreement. 
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APPENDIX A 


To simplify the solution of equation (16), an approximation is discussed 
in this section. 


The fol lowingrelation is given because (1 - F(t,a>) is a nonincreasing 
function against t. 


M 

n ' 

J-o 


y m 

1*2 


jl - F ^(m - 

3 >At,a f ) j 


* ( 

/ \) Y 


i-F 

f(3 - 3 )At 1 ,a f j 

(A-l) 

3=1 

TP 

• 

II 

3 


(A-2) 

t 

> 

r-H 

II 

3 I. 


(A-3) 


if , -r 


(A-4) 


Equation (A-l) is reduced to 
3 


, Y m 

V - F(! r* i f ) } Y n jl - F (<S - * II jl - F (<m - j)At,a f )| 

3=1 3=1 ’ 


13 

^ n jl - F ((13 - 3 >At r a f )j 
3 = 1 ( ’ 


(A-5) 


If 3 and y are chosen as {1 
following relation is given. 

(3 Y 

n jl - F ((3 - j>At 1 ,a f ) ( 

3 = 1 ' ’ 


F(t r ,a)} Y is nearly equal to one, the 



P(t r ) 


(A-6) 
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In this study, P(£ r ) is calculated by choosing 3 and 
F(t r ,a)} Y < 1 . 
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TABLE I. - LOG NORMAL DISTRIBUTION OF CRACK GROWTH RATE IN THE DIMENSIONAL MODEL 


[C = C c Mja(n)Bo g n+1 .3 



a = 4 C 

a < 

*c 


Median 

Upper confidence 
limit, 

90 percent 

Median 

Upper confidence 
1 imi t , 

90 percent 

(da/dt)y 

V^Ca 

2Ca 

V / 2Ca{a/a c ) 0 * 414 

2Ca c 

{da/dt ) l 

Ca/\/2 

Ca 

Ca c {a/a c )°- 4M /V5 

Ca 


TABLE II. - LOG NORMAL DISTRIBUTION OF CRACK GROWTH RATE IN THE 
NONDIMENSIONAL MODEL 


[C = Mjo(n) , 5 = a/0 av , a c = a c /D av .] 



a = a c 

a < 

*c 


Median 

Upper confidence 

Median 

Upper confidence 



limit, 


1 imi t , 



90 percent 


90 percent 

(da/dOy 

\/ZCZ 

2Ca 

V2Ca(a/a c ) 0 - 414 

2Ca c 

<da/dt) L 

Za/y/Z 

Ca 

Ca c (a/a c )°- 414 /V^ 

Ca 


TABLE III. - CONSTANTS FOR THE CALCULATION OF CREEP IN 304 
STAINLESS STEEL AT 923 K IN AIR 


Constant 

Val ue 

Note 

B 

1 . 37xl0~ 18 

e c = Bo n Norton's law 
e c : Creep strain rate [h“ 1 ] 

o: Stress [MPa] 

n 

7.1 

M J 

0.51 

For semi-circular surface crack 

a c 

0.1 mm 

— 

D av 

0.02 mm 

Average of grain boundary length 
between two adjacent triple points 

Dst 

0.005 mm 

Standard deviation of grain boundary 

length between two adjacent triple points 

C c 

0.0096 MPa -1 
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CRACK LENGTH 


LOG (a) 


(b) SCHEMATIC DIAGRAM SHOWING RELATIONSHIP BETWEEN CRACK 
GROWTH RATE AND HALF CRACK LENGTH. 


FIGURE 1. 





NONDIMENS I ONAL TIME, t 


FIGURE 3. - SCHEMATIC DIAGRAM SHOWING STOCHASTIC CRACK 
GROWTH. 
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CUMULATIVE PROBABILITY, F(t,a) CUMULATIVE PROBABILITY, G(a,t) 



NONDIMENSIONAL CRACK LENGTH, I 

FIGURE 5. - CUMULATIVE PROBABILITY OF CRACK LENGTH AT 
TIME t FOR A SINGLE CRACK, G(a,t). 



FIGURE 6, - CUMULATIVE PROBABILITY OF TIME WHEN A CRACK 
REACHES LENGTH a FOR A SINGLE CRACK, F(t,a). 












.3 .4 .5 

NOND I MENS I ONAL TIME, t = (C c a n >BO n n t 

FIGURE 8. - DETAIL DISTRIBUTION OF F(t,I> NEAR 


F(t,a) = 0. 
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CUMULATIVE PROBABILITY OF CREEP LIFE, P(t 



0 I 1 1 1 

.3 .5 .6 .7 

NONDIMENSIONAL CREEP LIFE, t r 


FIGURE 9. - CUMULATIVE PROBABILITY OF NONDIMENSIONAL 
LIFE PREDICTED BASED ON THE MODEL OF MICROSTRUCTUALLY 
SHORT CRACK GROWTH. 
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NOND I MENS I ONAL CREEP LIFE 


SURFACE AREA, S, mm 2 



FIGURE 10. - RELATIONSHIP BETWEEN NOND I MENS I ONAL CREEP LIFE AND NUMBER OF CRACKS. 



10 1 10 2 10 3 5x10 3 

TIME TO RUPTURE t r . h 

FIGURE 11. - COMPARISON OF CREEP LIFE BETWEEN PREDI- 


CTION AND EXPERIMENT. 
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og = 98 MPa 
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